Creating nuclear spin entanglement using an optical degree of freedom 
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Molecular nanostructures are promising building blocks for future quantum technologies, provided methods 
of harnessing their multiple degrees of freedom can be identified and implemented. Due to low decoherence 
rates nuclear spins are considered ideal candidates for storing quantum information while optical excitations 
can give rise to fast and controllable interactions for information processing. A recent paper (Physical Review 
Letters 104 200501) proposed a method for entangling two nuclear spins through their mutual coupling to a 
transient optically excited electron spin. Building on the same idea, we here present an extended and much 
more detailed theoretical framework, showing that this method is in fact applicable to a much wider class of 
molecular structures than previously discussed in the original proposal. 
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I. INTRODUCTION 

The controlled generation of entanglement is a crucial task 
in quantum computing, quantum cryptography, quantum er- 
ror correction and for other quantum technologies. Experi- 
mentally, controlled entangling operations have been demon- 
strated in a wide range of systems, such as for example pairs 
of photons 1 1 1, atoms [2 |, ions [3), as well as between an atom 
and a photon [4], an ion and a photon [5 1, among many others. 
However, a different approach typically needs to be used for 
each of these systems, so that a given control methods may 
not readily be transferred to a different physical system. 

Employing molecular systems as 'quantum hardware' J6l|2l 
offers the advantages of high reproducibility, chemically engi- 
neered system properties, and the potential for self-assembly 
into more complex functional units. The key challenge for 
nuclear spin-spin entanglement is achieving fast and switch- 
able control over the interactions between adjacent quantum 
bits (qubits). Several previous publications have proposed 
introducing a mediator spin, whose mutual coupling to the 
qubit spins provides a route for a controlled entangling op- 
eration ll8 UT2ll . In these studies the mediator spin is usually 
of the same type as the qubit spins and possesses spin 1/2. In 
contrast, many molecular structures possess optically-excited 
triplet states (i.e. states with spin 1 character) that could be 
used as mediators for nearby nuclear spins [13]. Motivated 
by recent experiments lTT4l IT31 we develop theoretical control 
strategies for this class of systems. 



II. MODEL 

For our model we consider a structure comprising three in- 
tegral components: two nuclear spin qubits labeled n and «' 
and a mediator spin system with an optical degree of free- 
dom. We denote the ground state of the mediator as |0) and 
the first excited state as \e). Further, |0) is assumed to be spin- 
silent, while \e) possesses an electronic spin 1 character. The 
spin 1 (quasi)particle is thus created upon optical excitation 
of the mediator system. Importantly, the two spin qubits do 



not directly interact with each other, yet both are coupled to 
the excited state of the mediator with an isotropic Heisenberg 
interaction as schematically depicted in Fig. [T] The Hamilto- 




ground state excited state 



FIG. 1. (Color online) (a) There is no coupling between the mediator 
and the two nuclear spins in the ground (or vacuum) state |0). (b) 
The excited state of the mediator (e.g. an electron-hole pair) has spin 
1 character and couples to both nuclear spins. 

nian for such a three particle system in an external magnetic 
field is then given by 

H = -COnS^ - ayS,,,/ + \e) ((a e S Zte + Oq) (e\ 

+ \e)(AS n -S e +A'S n ,-S e +DSl e )(e\ , (1) 

where denotes the Zeeman splitting of the two qubits 
and co e that of the central spin 1 particle; D is the uniaxial 
zero-field-splitting, A and A' are the isotropic Heisenberg cou- 
pling constants, and coo denotes the (typically optical) excita- 
tion energy of the mediator. Here S z n /„r i e and S„/„/ / e are the 
usual component and total vector Pauli spin operators respec- 
tively. Throughout this article we assume \D\, A <C CO e and 
(O e > 0. 

The Hamiltonian ([T| can describe many different nanos- 
tructures, in particular a range of optically active molecules, 
in which case the mediator ground (excited) state corresponds 
to the absence (presence) of an electron-hole pair. The transi- 
tion between these two states can be induced by a short laser 
pulse of frequency Oo; moreover, the excited state will de- 
cay naturally through spontaneous emission. By contrast, an 
NV~ -centre in diamond surrounded by two 13 C atoms lfT6l 
is an example of a system in which the spin 1 mediator is 
not susceptible to decay and is ever-present. Many of the re- 
sults discussed in this article are equally valid for this kind 
of system. However, the following discussion primarily fo- 
cuses on the example case of the molecular system consisting 
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of a functionalised Cgo molecule with two functional groups 
as depicted in Fig. [2] 




FIG. 2. (Color online) A functionalised buckyball with two satel- 
lite atoms (blue). In reality the functional groups are likely much 
larger than a single atom, but for simplicity we here only depict the 
two relevant atoms with a nuclear spin 1/2, for example 13 C, 15 N, or 
3 1 P. Note that the two satellites (or rather functional groups) are not 
necessarily both the same, and that the system may thus be 'asym- 
metric'. 

In what follows we will first we analyze Hamiltonian ([T} 
for two different cases: that of a symmetric and that of an 
asymmetric system. In order to understand the system prop- 
erties and dynamics, we will determine the eigensystem using 
perturbation theory. We next present suitable protocols for the 
controlled generation of nuclear spin entanglement for both 
cases. We conclude the paper with a brief summary and a 
discussion of our results. 

III. EIGENSPECTRUM AND EFFECTIVE HAMILTONIAN 

A. Symmetric system 

By definition the symmetric system consists of two nuclear 
spins with co n = (0„i and equal hyperfine coupling constants 
A = A'. In the presence of the electronic excitation, i.e. af- 
ter the application of a suitable laser pulse, the Hamiltonian 
governing the spin dynamics is given by 

H sym = (e\H\e) = ~(0 n S ZJ1 + (0 e S-. e ~ (O n S Zt „i 

+A(S B -S,+S B ,-S e )+DS| e , (2) 

where we have neglected the term (Oq that is proportional to 
the identity and unimportant for the dynamics. Since the elec- 
tronic Zeeman splitting co e is typically the largest energy scale 
of the system (and in particular much larger than a> n ), it is safe 
to assume that 

\(0„\ ,\D\ ,A co e . (3) 

Based on the above assumption, we employ degenerate per- 
turbation theory to determine the eigenspectrum of H sym . We 
begin by partitioning the Hamiltonian as follows 

^sym = ^O.sym "i" ^sym ! (^) 



where 

#0,sym = -OinS^n + (O e S zfi - CO n S z n i +DSj e (5) 

will be treated exactly, and the perturbation is given by 

H sym =A(S n -S e +S n '-S e ) . (6) 

AsHq S y m is diagonal in the computational basis (defined as 
{|eni« 2 > I e = T ± ,T and n l/2 =t4}) and Ut |#o, sy m| It) = 
(tl |#o.sym| t4-)> degenerate perturbation theory is required. 
However, since (TjY \H' sym \TiZ) = for j,l = — ,0,+ and 
Y,Z =|t>tl the degeneracy is not removed in the first or- 
der. To find the correct zeroth order wavefunctions, we need 
to solve the corresponding secular equations to second order 
ifTTl . Making use of the assumption ([3]) we then get the fol- 
lowing second order eigenenergies: 

E JA = I:,-. - V, E J%2 - -jco e + \j\D (7) 

Ej,3 = Ej t2 - Si E jt 4 = Ej.2 + Sj , (8) 

with the corresponding eigenvectors up to first order (see 
Fig. [3): 

\E ]A ) = \Tj U) \E j2 ) = -L (-\ Tj |t> + \Tj U>) (9) 
\E ]A ) = \Tj ft) \E Jt3 ) = i (\Tj It) + \Tj tl>) (10) 

and where we have defined the following quantities for com- 
pactness: 



e_ = co„ -A + 2a-, £o = «„ + 2a + , e+ = co„ +A , (11) 
8- = — 2a_, Sq = — 2«o, and 5+ = 2a + , (12) 




FIG. 3. Eigenspectrum of a symmetric system up to second order of 
perturbation theory (eigenvectors not normalized). Our basis states 
obey the following relations S~ e \T±) = =f|7±), S- f |7o) =0, S zn \ | 
} = [!>, and S z ,„[t>= -| t>- ' 

Based on this spectrum, we now approximate H sym in the 
following way to obtain an effective Hamiltonian: 

// sym «y t diag(£_ L1 ,...,£ 1 , 4 )y =:# sym ,eff , (14) 

where V is the matrix of the approximate eigenvectors (see 
Fig. [3} 

V = (|£_,i)|£_^)...|£ +(4 )) . (15) 
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In the computational basis H sym e g consists of three blocks 
with the following structure: 



V 



x x 

X X 



x denotes a nonzero entry 



(16) 



where each block belongs to one of the three spin states T , 
7b, or 7+. Correspondingly, we define the three subspaces 
3Tj spanned by the vectors 1 7) ||) , 1 7) |t> , I Tj U) , and | 7) ft) 
with j = — ,0,+. this paper, the index j will denote the elec- 
tronic spin state of the excitation whereas the letter i always 
indexes the four nuclear spin states. 

In writing Eq. ( fl4| i we have neglected all matrix elements 
between the subspaces and 2Tj (i ^ j); the effect of these 
matrix elements is negligible because the electronic Zeeman 
splitting co e is much larger than the nuclear Zeeman splittings. 
The dynamics in each of the 3~j is therefore closed and can 
be described by an effective 4x4 Hamiltonian of the form 
given above. We shall now analyze the nuclear spin dynamics 
arising from each of those effective subspace Hamiltonians. 
A term in the Hamiltonian that is proportional to the iden- 
tity has no effect on the dynamics within the subspace, so we 
can subtract D + co e + a_ , ao, and D — co e — a + from the sub- 
space Hamiltonians for 51 , and £f + respectively (thereby 
neglecting unimportant phases that are common to all states 
in each block). We can simplify the resulting Hamiltonians 
H sym e ff j further by transforming to a suitably chosen rotat- 
ing frame through the transformation 



dU) 



Hsym, eff, j — UjH sym , eff, jUj + iU j 



where 



and 



U j {t)=R ZA (<$> j t)®R z „,{tyjt) 



<j)± = — 0)„=fA — a± 
<j)Q = — CO„ — a_—a^ 
Rz,i{<t>) = expHS^) 



(17) 



(18) 



(19) 
(20) 
(21) 



In this rotating frame we obtain the following Hamiltonians 



^sym, eff, j 



(22) 



where k(j) = 1,1,-1 (with j = — ,0,+, as it will be through- 
out the paper). Each of these effective subspace Hamiltonians 
therefore induces the dynamics of a direct XY-coupling of the 
two nuclear spins, with a time evolution corresponding to Rabi 
flopping between the nuclear spin states | J,t) an d I tl)- Th e 
Rabi frequency aj depends on the spin state j of the optical 
excitation. According to Eq. ( fT3~| , we expect a + to be fairly 
similar to a_, whereas in comparison ciq will be much smaller. 
We will return to the symmetric case to discuss entangle- 



ment generation in Section IV A but first complete our analy- 
sis of the Hamiltonian for non-symmetric cases. 



B. Asymmetric system 

Next we consider an asymmetric system with unequal nu- 
clear Zeeman splittings and/or unequal nuclear - electronic ex- 
citation hyperfine coupling. The excited state Hamiltonian is 
then given by: 



77asym = (e\H\e) 



— CO n S z ,„ + (£> e S z e — (8 n tS z n i 

+AS n -S e +A'S n ,-S e +DS 2 zl 



(23) 



For the purpose of applying perturbation theory, we consider 
the nuclear-spin-mediator coupling as the perturbation and 
split the Hamiltonian as follows 

77o,asym = —(QnS z ,„ + (O e S z , e - (D n iS z rf +DS 2 ze and (24) 

^asym = ' + A S„/ • S e . (25) 

For the asymmetric system we shall assume that the param- 
eters A and A' and a>„ and ©,/ differ sufficiently to fulfil the 
following two inequalities 



1 



A' -A±{(D n t-(Dn)\ > |fl±|,|a±| 



(26) 
(27) 



where aj and a'j are as defined in the previous section, with 
the latter using (O n * and A' rather than co„ and A. Under these 
assumptions, non-degenerate perturbation theory can be used 
for calculating the eigensystem. The state corrections to first 
order perturbation theory then introduce a small mixing of 
computational basis states with corrections of magnitude 



1 



V2D±{co e - 



and 



1 



A' 



v / 2D±{(Oe + (O n i) 



(28) 



which are negligible for A, A' <§; co e . The eigenstates of 7/ asym 
thus coincide with those of 77o. asym (which are simply the com- 
putational basis states) to a very good approximation. 

Analogously to our approach in the symmetric case, we 
proceed by analyzing individual effective Hamiltonians for 
the subspaces 3Tj (j = — ,0,+). After adding a suitable con- 
stant to each subspace, we obtain three diagonal Hamiltonians 



-f^asym, eff,/' — tyjSz.n' + §jSz,n 



(29) 



where 0j, (j)j are as defined in Eqs. ( fT9] l and ( |2"0| ) and the prime 
denotes that (Q„> and A' are used for the expression instead of 
(0„ and A. All three effective Hamiltonians are local, meaning 
there is no direct coupling between the two nuclear spins. 



C. Crossover between a symmetric and an asymmetric system 

In a certain region of parameter space neither the symmetric 
nor the asymmetric analyses are justified. In this section we 
consider this 'crossover' regime, so that in the later discussion 
we will be able to interpolate between those two cases. 
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We start with the approximated symmetric Hamiltonian 
#sym,eff and treat the asymmetry as the perturbation H{. , i.e. 



^co — ^sym.eff "F ^co I 

H' QO = -AiCO„S z .„i +A 2 AS„/ -S^ 



(30) 
(31) 



where Ai is the fractional difference between the two nuclear 
Zeeman splittings, A; = — -, and A2 is the fractional dif- 
ference between the two coupling constants, A? = ^-j^ . 

We begin with the vectors \Ejj) and values Ej\ (j = 
— , 0, +; i — 1 , . . . , 4) from section III A as the eigenstates and 
eigenvalues of Hamiltonian // sym .eff ■ As in the other two cases, 
we can still neglect those terms of H co which couple different 
projections of the excitation spin with strength A2K/ \/2 since 
the relevant states differ in energy by about co e , yielding an 
effective crossover Hamiltonian i/ C o,eff with the eigenvectors: 



\E j , l ) = \T j U); \E h2 ) = ±=[a lll \T l \ [} + a rll \T, ]\|) 



(32) 



1^,4) - \Tj ft); \e j3 ) = (09,3,1 \Tj It) + awlTj tl>) ; 

(33) 



where 



a 



7.2,1/2 1 



a 



7.3,2/1 1 



1 ± sgn{aj)fj 
,a )'f] 



(34) 



and 

f± = l -(A 2 A±A x co n ) and f Q = -^Ai(0„ 
The corresponding eigenvalues are 

E jd/4 =E jA/4 ±k(j)fj 

Ej,2/3 = E j,2 + k(J) ( a j T sgn(a j)\Ja)+.ff} 



(35) 

(36) 
(37) 



where k( j) is again 1,1,-1 for j = — , 0, + respectively. 

It is easy to see that the eigenstates reduce to those of the 
symmetric case for a vanishing perturbation (fj = 0), whereas 
they tend to the computational basis states for a larger pertur- 
bation as required for the asymmetric system (|/)| S> \aj\)- 



A. Symmetric system: entangling time evolution 

For the symmetric system, we shall exploit the effective 
XY-coupling between the spin states | ft) and | for the 
generation of entanglement. The fact that the magnitude of 
this coupling depends on the spin state of the excitation will 
allow us to control the interaction. 

The free time evolution in any of the subspaces takes 
suitable initial product states of the nuclear spins to entangled 
states at certain subsequent points of time. In order to quantify 
the performance of the desired operation, we consider the en- 
tangling power of the unitary operator that describes the time 
evolution of the system. The entangling power is defined as 
the mean linear entropy produced by the unitary operator act- 
ing on a uniform distribution of all (pure) product states lfl8l . 
A maximally entangling two qubit gate, e.g. the CNOT- or the 
CPHASE-gate, possesses an entangling power of 2/9. 

For the symmetric system the entangling power of the free 
time evolution Uj(t) = exp(— iH' symsff jt) is given by the fol- 
lowing simple expression for each of the three subspaces 3~f. 



i(2ajt)) &in 2 (ajt) 



(38) 



The entangling power ej is only a function of the coupling 
strength aj between the states \Tj 4-f) and \Tj tl), becoming 
maximally entangling at times that are odd integer multiples 
of ?max = %l (2a j). However, the characteristic timescale over 
which entanglement builds up differs significantly between 
the subspaces, since a±/ao = co e /(2D) to leading order, and 
this ratio is typically much larger than one. Hence, <?o(?) ~ 
for t < %j (2a±) when |D| <C co e as we have assumed so far. 
Fig. |4] shows the entangling powers e + and eo for a typical 
ratio of a+/ao = 32. 
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FIG. 4. (Color online) Entangling power e +/ / of the time evolution 
operators [/ + /g(0- Here a + = 32aq (see text). 



IV. CONTROLLED GENERATION OF ENTANGLEMENT 

The contrasting analyses presented in the previous section 
suggest that the dynamics of the system will vary significantly 
depending on its parameters. Therefore, different approaches 
for achieving controlled entanglement between the two nu- 
clear spins will be required. In this section we will show how 
to do this in each case; we begin with the symmetric system 
introduced in section llH AI 



Based on these observations, we can formulate a protocol 
for the controlled generation of entanglement. Consider a sys- 
tem that is initially in its ground state (i.e. there is no excita- 
tion coupling to the nuclear spins). A laser pulse creates the 
excitation, which will in general be in a mixture of the three 
spin states \T-), \Tq), and \T + ). For the present discussion, we 
will assume it to be in the state | To) ; the discussion of an mixed 
initial state will follow later. A short microwave pulse allows 
us to flip the state of the excitation selectively to either the 
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|7+) or the |7L) state since the transition frequencies are split 
by the ZFS. The entangling dynamics proceeds much faster 
in these outer subspaces, reaching the first entangling power 
maximum after a time t = We can now apply another 
microwave pulse to flip the excitation back into the \Tq) state, 
where the further evolution is much slower and entanglement 
can be preserved. 

There are, however, a few subtle points worth pointing out. 
First, similar to the case of applying the maximally entangling 
CNOT operation, only suitable initial states will evolve into 
entangled states. Second, certain particular input states can 
reach maximal entanglement in less time than ^ ; for exam- 
ple, the time evolution U{-^-) takes \T + |t) to a maximally 
entangled state. 

Finally, the (optical) excitation, whose natural lifetime z 
is assumed to be longer than the time it takes to build up 
the entanglement ^ < T, needs to be destroyed quickly 
enough such that no further (slower) evolution unwinding the 
achieved entanglement occurs in the ,% subspace. This can 
either be achieved with a coherent optical ft -pulse, or alterna- 
tively by simply waiting for the system to decay back to its 
ground state if the following hierarchy of timescales exists in 
the system 



eT, 



11) -+ |0 ||), \eTj) j= (t| It) + 1 tl» |0) ^ (T| It) - 



% 71 

T~ <T<< tTTT 

2a± 2 flo 



(39) 



However, in contrast to the de-excitation with a coherent laser 
pulse, it is not immediately obvious that the nuclear spin en- 
tanglement can survive the optical decay process; we will 
therefore now take a small diversion to analyse this decay pro- 
cess in detail. 



Decoherence due to the optical decay process 

In general, the optical decay of the mediator induces deco- 
herence on the nuclear spins. In order to quantify this, we will 
make use of the quantum optical master equation (for a full 
derivation see e.g. Ref.fl9l: 

U(t) = £ e^'-^T{co)(A(co)p(t)A(co'y 

£0,6)' 



-A(co'yA(co)p(t)j +H.c. , (40) 

where p denotes the density matrix in the interaction picture, 
r(o') is the rate for a transition with frequency co and H.c. 
is the Hermitian conjugate. The sum is taken over all opti- 
cal transitions of the system with transition operators A(a>) as 
defined by 



A(a») = £ n(£)MI(£') 



(41) 



E'-E=w 



where E and E' are eigenvalues of H that differ by co and 
Tl(E) denotes the projection onto the eigenspace belonging to 
the eigenvalue E. £F denotes the system's optical dipole oper- 
ator. The symmetric system features twelve optical transitions 



| tl», and \eTj ||) -> |0 ||) for j = -,0,+. We make the 
additional assumption that the optical decay rates are all equal 
and can thus be characterised by a single optical lifetime z. 

Equation ( |40] ) can often be simplified by applying an in- 
stance of a rotating wave approximation (RWA), based on the 
assumption that fast oscillating terms average out [ 19], giving 
this more common form of the quantum optical master equa- 
tion: 



p(t) = L r (co){A(co)p(t)A(coy 



-A(o) t A(o)p(f)) +H.c. . (42) 



dt 



However, the RWA is only justified when \co — o'p 1 is small 
compared to the relaxation time of the system z. Under the as- 
sumptions of equation ([39]) this is fulfilled except for the two 
transition frequencies C0\ = £b,2 + (Oq and C0i = £0,2 + (Oo — 
Sq which corresponds to the transitions |e7o) (=f| It) + 

I tl» -> |0) (t| It) + I tl)) ■ Therefore we can safely ap- 
ply the RWA to all remaining frequencies, obtaining: 

~/>(0= E e'«°'-^nco)(A(co)p(t)A(co'y 
at w.m'es 

-A(a>'Ma>)p(t)) 
+ £ r(o)('A(a))p(f)A(o) + -A(a)) t A(a))p(f)) +H.c. 

(43) 

= £r(fl>) (A(a))p(r)A(to) t -A(co) f A(co)p(t)) (44) 

CO 

+ e 2 «°"^(2A(co l )p( t )A(co 2 y- 
A(fi) 2 ) t A((»i)p(f)-p(r)A(£i>2) t A(a)i)') +H.c. 
where S = {(Oi , ©2} and using phenomenological decay rates 



r(o>) 



^ if co > 
if co < 



(45) 



This definition of decay rates describes the typical situation of 
only spontaneous emission occurring, since stimulated emis- 
sion or absorption of photons are proportional to the negligible 
power density of thermally activated photons in the environ- 
mental modes. 

Rather than solving the above master equation for the entire 
Hilbert space of our system, we are here only interested in the 
4x4 density matrix p* of the two nuclear spins after the decay 
has occurred: 



(0|/J(f = oo)|0) 



(46) 



Limiting the following discussion to this nuclear spin sub- 
space, it is easy to see that only a small number of ele- 
ments of pf can be populated by the decay process. These 
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non-zero elements are p a h with a = b =||, a = b =ft> and 
a,b e {|t,tl}- The final nuclear spin state after the optical 
decay can be compactly written as 



PL = (eT a\po\eT b) + £ £ «rf)a|/Jo|«2) 



";=± 



<rf)a|/Jo|e7)i>) , (47) 



where po = p(0) is the full system's density matrix in the in- 
teraction picture before the decay process and the bar on top 
of the nuclear spin states denotes that these are flipped, i.e. 
| =t and j" =|. Hence the coherences between \eT± ||) and 
\eT± tl) do not necessarily survive, but all the coherences be- 
tween \eTo ||) an( J k^b tl) survive the optical decay, reflect- 
ing the fact that no RWA approximation has been made for 
the transitions |«r )(±| It) + I U» ->■ |0)(±| It) + I tl))- 
Physically, this means that these two transitions produce in- 
distinguishable photons due to overlapping emission spectra. 
In the absence of other decoherence processes the spectra are 
simply given by lifetime broadened Lorentzians |20|: 



L{co) 



1/T 



(C0-C0i) 2 + (1/T) 2 



fori = 1,2 



(48) 



The photons emitted in all other ten transitions can be in prin- 
ciple be distinguished, so that we are effectively dealing with 
eleven distinct, incoherent decay channels. It is worth not- 
ing that there are three decay channels which populate the 
|0)(±| ||) + I tl)) states, one each for the excitation spin 
projections 7o±, but only the decays from the ,% subspace 
preserve coherence. 

In Fig.|5]we plot the purity Tr(p(f) 2 ) = Tr(p(f) 2 ) to illus- 
trate the decoherence induced by the optical decay for three 
different initial density matrices p,(0) = |y/,)(i//, | with 

\Yi} = \eT + U}, \Y2) = \eT U) , and 

\ ¥3 ) = l - m (i id + 1 it) + 1 ti)+ 1 tt» • ,49) 



Solving the master equation ( |43| i analytically, we then obtain 
the following final density matrices of the nuclear spins after 
the decay: 



pf=kiitxiti+itixtii) 



(50) 



Pi = 2 + 2gg T 2 ( 8 ° A It) (It I + iSor\ It) (tl I (51) 



,/ -2 + 25 W"° l 

-idox\ t|)(|t| + (2 + 5 2 T 2 )|tl)(tl 
p{= + | |t)<tl I + ltl)Utl) , (52) 
with the corresponding purities 

Tr(pf) 2 = l/2 and Tr^) 2 = 3/8 , (53) 



Tr(p|) 



2^2 



f 2 _ 2 + g, 2 T 2 © 

' ' " 2 + 25 2 t2 ~ 2 



(54) 
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FIG. 5. (Color online) Evolution of the purity Tr(p(?) 2 ) of a system 
initialized in an excited state during the optical decay process. The 
three different initial states shown are defined in l |49| >, and zSq = 
0.05 < 71. 



At first it may seem surprising that an initial state \eT+) | tl) 
can end up in a complete mixture of |0) 1 1|) and |0) | tl) ■ This 
is due to the assumption ( |39l > which underpins the optical mas- 
ter equation ( |43| l. This means we have implicitly included the 
fast Rabi oscillations between | ||) and | tl) in the ST+ sub- 
space while the system is waiting for the decay. In contrast, in 
the ,% subspace the inherent dynamics is much slower so that 
the final result depends on the relative magnitudes of 5q and 

T. 

Perhaps surprisingly, decay due to the spontaneous emis- 
sion of a photon does not act as a source of decoherence if 
two conditions are met: (i) the system decays from the sub- 
space spanned by the two states |7q ||) and \To tl), and (ii) 
t8q <C it, i.e. the energetic splitting of these two states is 
small compared to the inverse natural lifetime. This property 
can be turned into a powerful feature for suitable molecular 
systems, which we shall exploit in the following. 



Dealing with a mixed electronic excitation 

So far, we have assumed that the creation process yields 
a completely polarized excitation, enabling the simple proto- 
col for the generation of entanglement described in a previ- 
ous section. Motivated by recent experimental data from a 
promising candidate molecule [ 14], we now analyze the im- 
plications of having an initial mixture of the states | To), and 
|7±). Experiments on a 13 C labeled methano-carbon of the di- 
ethyl malonate mono-adduct (DEMF) reveal that the popula- 
tion of the electronic excitation in a sample oriented along the 
z-axis are equally distributed between \T + ) and |7L). In this 
case the lifetime of the excitation also depends on the state of 
the excitation, being much shorter for the | To) state compared 
to|7±). 

In the following we shall demonstrate how the generic pro- 
tocol presented earlier can be adapted to accommodate for the 
properties of the specific system presented in Ref. 14 After a 
short laser pulse for the optical excitation, the subspaces 5fo,± 
were found to be populated as follows p_ = 0.49, po = 0.02, 
and p + — 0.49 with associated lifetimes T_ = 0.57 ms, To = 
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0.02 ms, and T+ = 0.57 ms. We shall assume that the nuclear 
spins are initialised in the state | J/|~) . 

The basic idea of putting the system into the \T±) states 
to let the free time evolution generate entanglement, followed 
by (mostly) switching off the interaction in the ,% subspace 
remains unchanged. As before switching between different 
electronic states is accomplished using microwave pulses that 
are fast on the timescale of the nuclear spin evolution. The 
adapted protocol proceeds in two stages following the optical 
excitation. First, we let the desired entanglement build up in 
the £? + subspace by waiting for the time t = We then 
swap the populations of |7q) and \T + ) and wait until the en- 
tangled populations have decayed. The difference in decay 
rates 1/fo 3> 1 ji± means that the population of |7L) largely 
survives once the majority of |7q) has decayed to the ground 
state. Second, population in the subspace is maximally 
entangled at times that are odd integer multiples of We 
pick the first such point of time after the \Tq) has emptied out 
and apply another microwave 7£-pulse to swap the populations 
of | 7q) and |7L). Once more, | To) will quickly drain into the 
ground state, meaning there is now no more excited popula- 
tion left. Ideally, we are left with an almost fully entangled 
nuclear spin state. 

However, the success of the above described protocol is 
predicated on the coupling strength A of the nuclear spins to 
the excitation. In particular, for a very small coupling strength 
A it takes a long time to entangle the nuclear spins t w ^f, 
so that there may be a substantial probability of the optical 
decay having occurred before the nuclear spins can become 
properly entangled. In this case the left-hand inequality of 
equation <j39j is violated. On the other hand, if the coupling 
strength A assumes very large values, then the right-hand side 
of the inequality ([39) is violated. In the latter case the pho- 
tons resulting from the transitions |e7b) 4= (=f| ]/[) + | f|)) — > 

|0)4|(=f| It) + 1 tl)) become distinguishable and the decay 
is hence no longer coherence preserving. In Fig. [6] we re- 
gard the hyperfine coupling A as a tunable parameter and plot 
the entanglement of formation OT1 of the final state of the 
two nuclear spins. We consider two different initial states for 
the mediator spin: a completely polarized state and the mixed 
state reported in Ref. Ifl4l . 
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FIG. 6. Entanglement of formation of the nuclear spins after apply- 
ing our protocols as described in the text for different initial polar- 
izations of the excitation: Ppoiarized = |7b){7()| (solid) and p crys tai = 
0.49|r_)(r_|+0.02|r )(r |+0.49|r + )(r + | (dashed). For the first 
case we used the simple protocol described first in this section; the 
switching time used is . In the second case we use the enhanced 
protocol described at the end of this section; switching times here 
are and . The nuclear spins for both curves are assumed to 
be initially in the state | 4-t); f° r the parameters we use the values 
found in a recent characterization experiment 1 14] D = —296 MHz, 
G) e = 9.6 GHz, Ohi = 3.7 MHz, T_ = 0.57 ms, To = 0.02 ms, and 
T+ = 0.57 ms. 



CO e ^> \(0„\. This allows us to assign a separate entangling 
power ej to the effective Hamiltonians describing each of the 
blocks. 

We now give an easily provable lemma which will enable 
us to write the relevant analytic expressions of the entangling 
powers ei for all three subspaces: 

Lemma. Let H be the time-independent Hamiltonian of two 
spin j particles with the following two properties: 

1. \E,) = | U),\E 2 ) = -a\ ;t> +b\ U),\E 3 ) = b\ It) + 
a\ tl), and \E 4 ) = | ft), with \a\ 2 + \b\ 2 = 1 and a, be K 
are eigenvectors ofH. 

2. The eigenenergies ofH satisfy E\ — Ei — £3 + £4 = 0, 
then the entangling power ofU(t) = exp(-iHt) is given by 



Robustness to imperfections in the symmetry of the system 



,/3) = ^(^-^sin 2 (§)-^(^-^) 2 sin 4 ( 



(55) 



The previously described protocol for the controlled gener- 
ation of entanglement assumes a perfectly symmetric system. 
In the following we will analyze the degree of imperfection 
in the symmetry that may be tolerated. We already know that 
the crossover between the symmetric and the asymmetric case 
is not entirely abrupt. In Section [HIC| we have found expres- 
sions for the eigenstates and eigenvalues that can fully inter- 
polate between the symmetric and the asymmetric case. In 
this crossover case, the effective Hamiltonian still consists of 
three distinct blocks each corresponding to the spin state of the 
excitation \Tj) (j = — ,0,+), and matrix elements connecting 
the blocks are negligible due to the large energy difference 



where /3 = \E-$ —E2\t. In addition we have e(a,j3) = e(£>,/3). 

Applying the above Lemma to our effective Hamiltonians 
for the crossover case gives three expressions e) = e(cj,Pj) 
with 



fj 



Pj=2tJaj+fj and cj=- 1 



(56) 

It is easy to see that the parameter aj now directly competes 
with the strength of the perturbation fj in the expression for 
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the entangling power. With the following measure for the 
asymmetry 



Xj 



(57) 



we can write the entangling power compactly as 



3 + 4j# + cos [2ajtJl + xj ))sin[ajtJl + %] 



9 \}+Xj 



(58) 



We plot ej for different asymmetries % in Fig. [7] Not surpris- 




FIG. 7. (Color online) Entangling power ey of the free time evolution 
for an asymmetric system, when the excitation is in the state \ Tj) for 
different strength of the asymmetry %j. 

ingly, the entangling power decreases with increasing asym- 
metry Xj- The maximum of the entangling power, which is 
achieved for /3 ; = kn where k is an odd integer, is given by 



9 ( c ; c j> 9 \ c j c j) 9(l+ X ]) 2 



; (59) 



this expression is plotted in Fig. [8] There is no significant 




FIG. 8. (Color online) Maximally attainable entangling power nij of 
the free time evolution, when the excitation is in the state \Tj) with 
respect to the strength of the asymmetry Xj- 

reduction in the achievable entanglement power as long as 
Xj < 1/2, but the maximum drops quickly outside this regime. 
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Ax (10- 4 ) 

FIG. 9. Maximal attainable entangling power m_,mo, and m + as a 
function of the asymmetry in the Zeeman splitting A] = ° > "' m " and 



hyperfine coupling asymmetry A2 = A A A . Other parameters are: 
A = 2.5 MHz, co„ = 3.7 MHz, D = -296 MHz, co e = 9.6 GHz. The 
entangling power mj only reaches its maximum limit of 2/9 for cer- 
tain values of Aj and Aj. See main text for a more detailed discus- 
sion. 



We note that for equal hyperfine couplings (A2 
equal nuclear gyromagnetic ratios 



X± 



a± 

N 



03 e 

2\b\ 



> 1 



0) but un- 



(60) 



meaning that the ,% subspace's entangling power eg 1S much 
more affected by the asymmetry than e±. Fortunately, the dy- 
namics in this subspace is also the slowest, so that it can still 
conveniently serve as a shelf for entangled states generated 
in or ^1 until the optical excitation has been de-excited 
or decayed. Importantly, Eq. |59| implies that our scheme 
is robust against the small deviations from a perfectly sym- 
metrical system which one might expect in real-world experi- 
ments. Further, intentionally introduced small differences be- 
tween the frequencies co n and co n i (e.g. a chemical shift caused 
by different surrounding environments) may actually be use- 
ful for individual control and tomography of the nuclear spins, 
while a high-fidelity entangling operation is still possible. 

So far we have discussed the behaviour of the entangling 
power in terms of the parameter Xj- While the dependence of 
the nij on Xj is universal across the three subspaces we ob- 
tain qualitatively different results when considering plots that 
are based directly on the A1/2 asymmetry parameters. Fig. [9] 
shows the entangling power as a function of Ai and A2. We see 
that the behaviour is indeed quite different in each of the three 
subspaces: In the subspace we obtain a ridge along which 
an asymmetry between the nuclear Zeeman splittings and the 
hyperfine coupling constants completely cancels out, meaning 
a perfect operation is possible even for a system that is quite 
far removed from being symmetric. In contrast, the asymme- 
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tries add up in the subspace, so that the error tolerance is 
much reduced in this case. Finally, the ,% subspace is only 
sensitive to Ai, i.e. the difference in the nuclear Zeeman split- 
tings without any dependence on the hyperfine constants (see 

The entangling power vanishes completely in the limit of 
an entirely asymmetric system (which we take to be defined 
by the inequalities |26) and (p7|)). In this case %j^> 1, and the 
eigenstates consequently coincide with the computational ba- 
sis states, and the eigenvalues are such that the free time evo- 
lution no longer generates any entanglement (i.e. requirement 
2 of the Lemma is satisfied). We shall analyze this situation 
in the following section. 



B. Control methods for the asymmetric system 

For an asymmetric system we cannot rely on the system's 
free time evolution for the generation of entanglement; this is 
a direct consequence of the Hamiltonian being decomposable 
into local Hamiltonians (see Eq. (|29]i). Therefore, we need to 
apply a suitable sequence of radio-frequency and microwave 
control pulses to accomplish our aim of creating an entangled 
nuclear spin state. Hence, we proceed by analysing the dipole- 



allowed transitions of the asymmetric system (see Fig. 10 1. 

The asymmetric system possesses six (different) nuclear 
spin transitions on the radio-frequency scale, one per nuclear 
spin for each of the three spin states of the excitation. Re- 



ferring back to section III B we obtain these from the second 
order eigenenergies: 

o\i.j = | -jA -(>)„- a- (Sj- + 8j,o) - a+ (<5/,o + 8 jt+ ) \ , 

(61) 

Vrfj = | -J A ' - <°n> ~ a - i 8 j- + 8 j,o) - a' + {8j.o + 8j t+ ) | , 

(62) 

where S^j is the Kronecker delta, and as before j indexes the 
spin state of the excitation |7)) . Further, £% ; denotes the tran- 
sition frequency of the first nuclear spin and ft) r ' f ■ the transi- 
tion frequency of the second nuclear spin. In general all six of 
these frequencies may be distinct. 

Conversely, the spin state of the electronic excitation can 
be flipped conditional on the nuclear spin state using a suit- 
able microwave pulse. With four nuclear spin states and 
two excitation spin transitions, this gives a total of eight mi- 
crowave frequencies taking the excitation from \Tj) to \Ty) 
with (j,f) = (+,0) or (j,f) = (0,—) or vice versa. These 
are: 



CO 



CO 



Mw,U/tt 

W 
>w4t/tl 



CO e + {-\yD±-(A+A') 
tO e + {-\yD± l -{A-A r ) 



(63) 



(64) 



Here we have neglected second order perturbation theory 
shifts proportional to aj and a'-, as these are typically very 
small when compared to A and D. 

Several possibilities exist for exploiting this rich transition 
spectrum in order to realise an entangling operation. In the 




FIG. 10. (Color online) Microwave (blue) and radio-frequency tran- 
sitions (red) of the asymmetric system. The computational basis 
states are eigenstates of the system Hamiltonian (assuming Eqs. \26\ 
and \27\ are fulfilled). 



following we shall discuss three methods in more detail: a 
single microwave 2^-pulse, a pulse sequence of radio- and 
microwave pulses, and finally an adiabatic following method. 



1. CPHASE-gate through a selective 7.11-pulse 

Simply applying a selective 271-pulse with the frequency of 
any of the microwave transitions given in Eq. ( |63"j ) and ( |64| ) 
naturally implements a CPHASE-gate by imparting a phase 
of e m = — 1 to only one of the four nuclear spin states l22l . 

If the lifetime of the transition was infinite (and in the 
absence of other spin dephasing mechanisms), the 2^-pulse 
could be made perfectly selective, achieved by a pulse that is 
long in the time domain and accordingly spectrally narrow in 
the frequency domain 11231 . However, in practice the optical 
lifetime will be finite and this may limit the selectivity and 
thus the amount of entanglement that can be achieved. The 
challenge is to find the right balance between a fast pulse that 
is only partially selective, and a slow pulse during which the 
system suffers from the decoherence induced by the decay. 
We proceed by analysing the trade-off that arises from these 
constraints in the following. 

Suppose we are given an initial state that is an equal super- 
position of the computational basis states and a fully polarized 
state of the excitation, 



I initial) = \\To)( III) - 



lt) + ltl> + ltt» 



(65) 



and then apply a 2^-pulse microwave pulse with power Q.q 
and with a frequency cod corresponding to the energy dif- 
ference between the levels \T + "ft) an d \To tt)- To describe 
the dynamics of the excited system we use either the effective 
asymmetric or the more general effective crossover Hamilto- 
nian, whichever is more appropriate for the precise combina- 
tion of system parameters in question. In particular, //asym.effj 
is adequate whenever the eigenvectors closely coincide with 
the computational basis states, whereas H co & ff j is used oth- 
erwise. We apply the following criterion for discriminating 



10 



between the two cases: 



Heff.j 



J ^asym,eff J 

I #co,effJ 



|«/,2,i| 2 < 0.001 or 1^,2,1 1 2 > 0.999 
0.001 < \a jX \ 1 2 < 0.999 , 



(66) 

with Of/,2,1 as defined in Eq. ( [34| i. On top of the effective sys- 
tem Hamiltonian H e R we also need to model the microwave 
pulse (in the usual RWA), so that the total Hamiltonian during 
the pulse is given by: 



Huw = \0)(-(OnS z , n - 



-(O nl S z y)(0\ + \e)H ef{ (e\ 
+ \e) (£l Q S x . e - a D S z e )(e 



(67) 



As this Hamiltonian is time-independent we can use a quan- 
tum optical master equation like the one defined in Eq. ( |4"3| l 
to model the decay of the excitation in the interaction picture. 
The transition operators A (to) are as defined by Eq. ( pT} , with 
appropriate projectors onto the eigenspaces of H^. In our 
calculations we perform the RWA in the master equation (re- 
member that this RWA is different and independent from the 
RWA for the driving) whenever the difference of two frequen- 
cies differs by more than 30t _1 , where T denotes the lifetime 
of the excitation. 

While optical decay during the application of the pulse may 
preserve some coherence between the nuclear spin states for 
specific parameter combinations, this is no longer the case 
once the pulse has finished: the decay to the ground state 
for an asymmetric system in the absence of microwave driv- 
ing invariably destroys the nuclear coherences. Therefore, we 
must use a different approach for taking the system back to 
its ground state. The first possibility is de-excitation using a 
resonant optical 7C-pulse. Another option would be to signifi- 
cantly 'speed up' the decay process, e.g. by exciting the sys- 
tem into a different metastable excited state which is known 
to quickly decay to the ground state. Provided the lifetime 
of this metastable state is short enough, the wave function 
of the emitted photon does not carry information about the 
nuclear spin state, so that the nuclear spin coherence will be 
preserved. We shall assume that such a coherence preserving 
de-excitation can be accomplished. 

As mentioned above, the system is susceptible to poten- 
tially harmful decay events during the application of the 2n- 
pulse, such that the final nuclear spin state p nuc is a mixture 
of population that has spontaneously decayed and the remain- 
ing population to which the control sequence has been fully 
applied. Since we are now dealing with mixed states p nuc , 
the entangling power is no longer a suitable measure for the 
quality of our operation, and as before we employ the entan- 
glement of formation ll2D as an alternative benchmark. 

Assuming a simple top hat pulse profile in the time domain, 
the pulse duration for a 27T-pulse is t = where Ho is the 

applied microwave power. For optimal performance the right 
balance must be found between pulse selectivity and duration 
for each combination of system parameters and lifetime T. We 
thus maximize the achievable entanglement of formation by 
varying ilo to obtain 



EF* = maxEF(p n 
n 



As an example we choose T = 10 |xs to plot the quantity EF* 
in Fig. 11 as a function of A and A'. Larger hyperfine coupling 



constants allow a faster selective pulse, and the optimized en- 
tanglement of formation of p nuc hence increases with A and 
A'. Remarkably, even for a lifetime as short as 10 |os, a high 
entanglement of formation can be obtained with only moder- 
ate hyperfine coupling strengths. 

Finally, we note that the currently presented protocol also 
works for the symmetric system, where the entanglement op- 



eration then only takes a time t = sa ^Sp, which is much 
faster than our protocol discussed in Section IV A For a short 
optical lifetime this approach may thus be advantageous as- 
suming a spectrally narrow highly selective 2k pulse can be 
implemented. 
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FIG. 11. (Color online) (a) Maximized entanglement of formation 
of he nuclear spins after the 2^-pulse with respect to the two hyper- 

of 



111 



fine coupling strengths A and A', (b) Optimal duration t* 

the microwave 2^-pulse with respect to the two hyperfine coupling 
strengths A and A'. We use the illustrative lifetime X = 10 ixs while 
the other parameters used in this plot are motivated by Ref. 1151 : 
D = -320 MHz, co e = 9.7 GHz, a„ = 5.97 MHz, a n , = 14.74 MHz. 



(68) 
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2. Combined microwave and radio-frequency pulse sequence 



the ZFS were measured as 



In the previous section we have discussed a simple imple- 
mentation of the CPHASE-gate that is maximally entangling 
for suitable system parameters. However, other methods for 
creating maximally entangled states also exist, and these may 
be more advantageous if a particular final state is required. 
Here we shall briefly discuss an alternative control method 
that employs a sequence of microwave and radio frequency 
pulses instead of a single microwave pulse. Let us assume we 
have the initial spin state 1 7b ||) an d would like to create the 
entangled Bell state 4j(|7b ||) +\T tt)). It is impossible 
to perform this operation with only radio-frequency pulses. 
However, we can achieve our aim by 'shelving' parts of the 
population in one of the other electronic spin states. An exam- 
ple of how this approach works in detail is depicted in Fig. 12 
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A«6MHz, A' « 11 MHz, and D rj —320 MHz (70) 

at an external magnetic field of B = 0.346 T. Setting an upper 
bound of 0.01 for the unwanted transition probabilities, the 
four required pulses can all be applied in less than a 1 (xs for 
typical nuclear gyromagnetic ratios. The duration of the pulse 
sequence is thus short compared to expected optical lifetimes 
in candidate molecules, so that a high fidelity entangling op- 
erations using this protocol should be feasible. However, we 
refrain from performing a more detailed analysis as in the pre- 
vious section since the results are similar and little additional 
insight is gained. 



3. Implementation of a CPHASE operation with adiabatic 
following 

The last method discussed in this paper for creating en- 
tanglement in our system relies on the adiabatic following 
of system eigenstates, similar to the protocol described in 
Refs. IPT21 |24l . Here, it is implemented by slowly modulat- 
ing the intensity of a microwave pulse that is close to res- 
onance with one or several of the microwave transitions of 
the excitation spin. Prior to the application of the pulse, the 
(asymmetric) system is prepared to be in a superposition of 
computational basis states as follows: 

I initial) = I To) (fl 1 1 ||) + fl 2 1 It) + «3 I U) + «4 1 tt) ) , (71) 



FIG. 12. (Color online) Pulse sequence for the creation of the Bell with normalisation Yy 1 



state ^ (| To ||) + \Tq tt)) ■ The black energy levels are populated 
and the grey levels are empty. The dotted arrows denote a coherence 
between the adjacent energy levels. Only relevant spin levels are 
shown. 



Of course, the microwave and radio-frequency pulses also 
need to be sufficiently selective for this approach, which im- 
poses a minimal overall pulse sequence duration. Once more, 
we refrain from discussing sophisticated pulse shaping tech- 
niques [23 1, instead considering considering Rabi's instruc- 
tive formula for the transition probability of a driven two level 
system 



1 . Starting from this state, 



(0 



Q 2 



£2 2 +A 2 ' 



V£i 2 + A 2 r) 



(69) 



where £2 denotes the strength of the pulse and A the detun- 
ing. We apply this to the closest neighbouring transitions 
of the (resonantly driven) target transition and demand that 
£2 2 /(£2 2 +A 2 ) < 1 for all neighbouring transitions, allowing 
a rough estimation of the duration required for each pulse in 
the sequence. 

We take the example of a recent characterization experi- 
ment on an asymmetric system (phosphine oxide fullerene 
DMFPH) 03], where the two nuclear spins of interest are pro- 
vided by a hydrogen and phosphorus atom in the functional 
group attached to the fullerene. The hyperfine coupling and 



the microwave power is then varied such that adiabatic fol- 
lowing of instantaneous eigenstates occurs. Once the power 
is decreased again all population returns back to the compu- 
tational basis. During the pulse, the eigenstates are energet- 
ically shifted and thus pick up a dynamic phase. However, 
the precise shifts of the states differ due to the the hyperfine 
coupling, so that in general each of the four states acquire a 
different dynamic phase. This gives rise to an overall combi- 
nation that can be non-trivial and entangling l24l . 

Consider applying a microwave field with frequency cod 
whose power envelope is changed gradually following a Gaus- 
sian function £2(f) = £2oexp[— (f/ff) 2 ]. We apply this Gaus- 
sian microwave pulse from an initial time t = — 3<7 to t = 3<T 
and choose the frequency C0d in such a way that the pulse is 
off -resonant with all microwave transition in the system. The 
diagonal form of the system Hamiltonian // aS y m ,eff permits a 
description of the dynamics of each state in the superposition 
of Eq. ( |7T| with a Hamiltonian connecting all three excitation 
states of the form 





£l(r) 
V2 


\ 


a(t) 





Sl(t) 




v/2 


v° 


£2(0 

V2 


A i,2 J 



(72) 



written in the basis {|T_i), \T i), \T + i)} for / = {||,|t,tl,tt} 



with detunings An = ft> 



0o- 



COo and A, 



ft) 



+o0 



(0 D . 
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Note that the usual RWA has been performed. Choosing 0)d 
such that for each i 



a(0 



< A 2 ,i 



(73) 



then enables us to write an approximate Hamiltonian for each 
of the nuclear spin states (valid to first order in perturbation 
theory) as 




(74) 



To achieve adiabatic following the eigenenergies need to be 
varied slowly to suppress Landau Zener transitions between 
different eigenstates. Following Ref. [T2l this can be accom- 
plished under the following conditions: 



Go 
A? 



<d for / = {||, |t, tUtt} 



(75) 



As we have mentioned above, the net effect achieved by the 
adiabatic pulse is the acquisition of a phase 0, for each of the 
nuclear spin states. The 'right' combination of control param- 
eters <7,i2o, and A\j gives rise to an operation that is locally 
equivalent to a CPHASE-gate if the following condition is met 
G3: 



01 - 02 - 03 + 04 



(76) 



In terms of the evolution of the nuclear spins, the state | ^initial) 
has then evolved to 



\T Q )( ai e^ 1 1|) +a2 e^\tf)+a^\n)+a*e w <\ ft)), 

(77) 

corresponding to a local phase on each nuclear spin in addition 
to the application of a CPHASE-gate. 

We shall now address the question of how the optimal com- 
bination of control parameters may be found. The dynamical 
phase that is acquired by a state during the pulse is directly 
determined by the eigenenergy of the state \Tq{), yielding 




Imposing the condition ( |76| > and solving for a 
yields 



(78) 

2exp(-2x 2 )c/x . (79) 
ct(<Od,Q ) 




(80) 



where the sum is taken over the four nuclear spin eigenstates. 
To mitigate the effect of decoherence caused by the decay of 
the excitation we minimize the duration of the pulse under the 
constraint that the conditions in Eqs. |73|, f75] l, and |76) are 
fulfilled, thus obtaining an optimal a*. 



We incorporate the decoherence caused by the optical 
decay by modelling the time evolution with the following 
Schrodinger picture master equation 



dt 



p(t) = -i[J% op (o*),p(t 



+ £r(co) (2A(a))p(f)A(to) t - {A(fl>) + A(a>),/>(f)}) (81) 

in Lindblad form |[T9l . r(co) is as defined by E q. |45| ) and 
the Lindblad operators A((o) are determined by (|4l|i, where 
the projectors project onto the eigenspaces of // a sym,eff instead 
of J^pp- This simplification gives twelve (constant) incoher- 
ent decay channels, rather than considering time-dependent 
Lindblad operators and re-evaluating the validity of the RWA 
at every instance in time (as would be required for the time- 
dependent Hamiltonian). Effectively, our approach then over- 
estimates the destructiveness of the optical decay, thus giving 
a lower bound for the entanglement of formation of the nu- 
clear spins. 



Fig. 13 shows the results of a simulation that applies such 
an optimised Gaussian pulse with a pulse duration a* . For 
weak coupling strengths A and A', this approach achieves a 
somewhat lower value of the entanglement of formation than 
the dynamic 27T-pulse discussed earlier. However, a similarly 
high fidelity entangling operation is possible for stronger hy- 
perfine coupling. The term 'adiabatic following' can invoke 
the impression that the desired operation will be much slower 
than a dynamical implementation. It is therefore astounding 
that our adiabatic pulse only takes about twice as along as the 
dynamic 27C -pulse. Finally we note that the adiabatic method 
(where the pulse is applied off-resonantly rather than having 
to hit a specific resonance) is inherently robust against pulse 
imperfections. This could be a significant advantage for ex- 
periments with ensembles of identical molecules. In this case, 
static and driving field inhomogeneities will inevitably lead 
to an over- or under-rotation of some of the ensemble spins 
when a dynamical pulse is applied (leaving the wrong exci- 
tation subspace populated), whereas the adiabatic approach 
ensures that all populations end up in back in the correct spin 
levels. 



V. SUMMARY 

In this paper, we have given a detailed analysis showing 
how a transient optically excited state can be harnessed for the 
controlled generation of entanglement between two remote 
nuclear spins. We have identified control methods applicable 
over a wide range of system parameters and studied their per- 
formance with regard to the predominant decoherence mech- 
anism. For the symmetric system consisting of two identical 
nuclear spins as qubits, the free time evolution is naturally en- 
tangling, but the characteristic timescale of the dynamics de- 
pends on the state of the excitation and can vary over several 
orders of magnitude, opening up the possibility of effectively 
switching the interaction on and off. For an asymmetric sys- 
tem, a different route needs to be taken and we have presented 
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FIG. 13. (Color online) (a) (b) Optimal duration of the adiabatic 
pulse 6(7* optimised over a)£> and Q.q. (b) Entanglement of for- 
mation after applying a Gaussian pulse whose duration is charac- 
terised by ct*. As in Fig. 1 1 1 1 other parameters are: X = 10 [is,D = 
-320 MHz, co e = 9.7 GHz, fi)„ = 5.97 MHz, m n , = 14.74 MHz. 



one adiabatic and two dynamic methods for creating entan- 
glement in this case. We have also included a discussion of 
the crossover regime between the asymmetric and symmetric 
system to establish the robustness of the symmetric operation. 

We have shown that the symmetric control method can be 
remarkably robust against uncertainty or fluctuations in the 
coupling constants and nuclear Zeeman splittings. As another 
advantage of the symmetric system, the system can decay 
back to the ground state without destroying the nuclear spin 
coherence. Conversely, for the asymmetric system additional 
control is required for the de-excitation step, yet it is easier to 
address the nuclear spins individually for single qubit opera- 
tions, initialisation and read-out. 

Interestingly, the active control methods proposed for the 
asymmetric system are much faster than waiting for the free 
time evolution in the symmetric case, and they can also be ap- 
plied to the symmetric system if a short optical lifetime makes 
this approach advantageous. Finally, we note that the adia- 
batic control method is intrinsically more robust against con- 
trol pulse and static field inhomogeneities, making it uniquely 
suitable for experiments with ensembles of identical systems. 
Astonishingly, the time required for such an adiabatic opera- 
tion is only about twice as long as for its dynamical counter- 
part. 
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